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Abstract 

Thermal conductivity of graphene nanoribbons (GNR) with length 106 A and width 4.92 A after 
isotopic doping is investigated by molecular dynamics with quantum correction. Two interesting 
phenomena are found: (1) isotopic doping reduces thermal conductivity effectively in low doping 
region, and the reduction slows down in high doping region; (2) thermal conductivity increases 
with increasing temperature in both pure and doped GNR; but the increasing behavior is much 
more slowly in the doped GNR than that in pure ones. Further studies reveal that the physics of 
these two phenomena is related to the localized phonon modes, whose number increases quickly 
(slowly) with increasing isotopic doping in low (high) isotopic doping region. 
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I. INTRODUCTION 

The thermal conductivity is a interesting property of graphene. Recent experiments have 
measured the thermal conductivity to be an extremely high value. [lJjJJLThermal conducti vity 



can be theoretically studied by quantum molecular dynamics (MD)[3|, |4j, classical MD|5Hll|. 



14, classical MD [5 



or Klemens perturbation theory approach [121]. For the classical MD, the temperature should 

be corrected to account for the quantum effect on the thermal conductivity. This quantum 

r~i 

correction has been shown to be important below room temperature for nanotube[9( and 



nanowires lfj . Meanwhile, it has been demonstrated that the umklapp, defects and edges 
scattering have important effect on the thermal conductance in single layer graphene. [13] 
Isotopic doping also has very important effect on thermal transport. For nanotube, both 
experimental jl4| and theoretical J8[ results show that the isotopic doping can reduce the 
thermal conductivity remarkably. 

In this paper we investigate the isotopic doping effect on the thermal conductivity in 
graphene nanoribbons (GNR). The thermal conductivity of GNR is studied by classical MD 
with quantum correction. Our studies display two phenomena: (1) in low isotopic doping 
region, the thermal conductivity decreases rapidly with increasing doping, and 10% doping 
can yield 50% reduction of the value of thermal conductivity; while in the high isotopic 
doping region, the thermal conductivity decreases slowly with further increasing doping. (2) 
thermal conductivity of both pure GNR and isotopic doped GNR increases with increasing 
temperature, but the increasing behavior in doped GNR is slower than that in pure ones. 

The underlying mechanism for these phenomena is the localization of phonon modes. (1) 
For the first phenomenon, we find that a single isotopic doping center can localize three 
phonon modes with the frequency as 634.4, 1373.3 and 1629.9 cm -1 . Due to their localizing 
property, these modes are unfavorable for thermal transport. So generally the thermal con- 
ductivity will decrease after isotopic doping. We further show that the number of localized 
modes increases linearly with increasing doping in the low doping region, while this number 
increases very slowly with increasing doping in the high doping region. This mechanism 
results in the first phenomenon. (2) For the second phenomenon, at higher temperature 
more phonon modes will be excited and the quantum correction is less important. These 
two positive factors facilitate thermal transport. So generally the thermal conductivity will 
increase with increasing temperature. However, in doped GNR, more localized modes will 



also be excited at higher temperature, leading to negative effect on thermal transport. This 
negative factor competes with the former two positive factors, slowing down the increasing 
behavior in isotopic doping GNR. It should be noted that the increase of the thermal con- 
ductivity for graphene nanoribbon is true for nanoribbons where the transport is dominated 
by the edge effects. However, in large flakes of graphene (larger than the phonon mean free 



path), the Umk 



temperature. [13| 



app scattering will lead to decreasing thermal conductivity with increasing 



The rest of this paper is organized as follows. In Sec. II, thermal conductivity of graphene 
is calculated by MD, and results are interpreted by localization theory. Sec. Ill is the 
conclusion. 

II. MOLECULAR DYNAMICS SIMULATION 
A. simulation details 

The thermal conductivity is calculated through the direct MD method with the second 



generation Brenner inter-atomic potential applied|15|]. We apply periodic boundary con- 
dition in the width direction and fixed boundary condition in the length direction. The 
temperature for the left and ri ght part is controlled to be constant value Tl and Tr by the 



Nose-Hoover (NH) heat baths. 16 



171 ] The dynamical equations for the NH heat bath are: 



dn Pi_ 
dt vfii 
dpi dV 

*- = -*-- "'"" (1) 

Q = gk B Tr 2 , 

where rji = if atom % is not in the heat bath, and g is the number of degrees of freedom 
for atoms in the heat bath, r is the relaxation time of the heat baths and kept as 0.4 ps in 
our calculation. From the above dynamical equations, we can derive the exchanged energy 
between the heat bath and the system: 

pto+At 2 

E nh = - J2^dt, (2) 

J to i m i 



where the summation of % is taken over atoms in the heat bath. So the heat current is 181 ]: 



J = -^. (3) 

s At w 

In the steady state, the magnitude of heat current from the left and right heat baths should 
be equal: Jl = —Jr- So we can use J = (Jl — Jr)/2 to calculate the heat current of the 
GNR. The difference between left and right heat current, i.e dJ = \Jl\ — \Jr\, is used to 
determine whether the system has achieved steady state or not. Due to the numerical error, 
the value of dJ in the steady state is not strictly zero. This nonzero value is used to estimate 
the error bar for the value of thermal conductivity results throughout this paper. 

We use 8.5 ns for the system to relax and do average with another 8.5 ns. Once the 
temperature gradient is obtained from this temperature profile, the thermal conductivity is 
then calculated by the Fourier law: 

J = " K d? (4) 

B. quantum correction 



Since the Debye temperature Td of graphene is very high (above 1000 K).[19l42l| it is 
necessary to consider the quantum effect for the thermal conductivity. We do quantum 
correction for the results of the classical MD in the following ways. Firstly, the quantum 
temperature T q is calculated by equaling the total ensemble energy to half of the total 
phonon energy: [5| 

3 1 f°° 

-Nk B T md = - g(u)n(uj,T q )}Twdu, (5) 

* ^ Jo 

where T mi i is the MD simulation temperature. g{u) is the density of states (DOS). We 

calculate g(oS) from the Brenner empirical potential implemented in the 'General Utility 



Lattice Program' (GULP)[22|. n(oj,T q ) is the distribution function: 



nM) = e X p(hu/l B T q )-l (6) 

The thermal conductivity with quantum correction can be obtained from the classical MD 
results through following relationship [5j 

K>q\-Lq) = K md(,-'mdJ ry-p i (J) 



where K m d is the MD result at temperature T m( j and n q is result with quantum correction. 
The correction factor ( ™ d ) is shown in Fig. [TJ It shows that in T m d < 300K region, this 
quantum correction is more important, while the correction is small in the high tempera- 
ture region. This tendency is qualitatively consistent with the previous results on carbon 
nanotubes[9] and silicon nanowiresfll]. 

C. results and discussion 

At T m d=300 K, the value of the thermal conductivity of pure GNR in our calculation is 
about 61 W/(mK). The GNR has 400 carbon atoms, with length 106 A and width 4.92 A. 
This value is comparable with that of the carbon nanotube with similar length 9]. These 
small absolute values of the thermal conductivity obtained in theoretical works are the results 
of the very small size of the examined graphene nanoribbon compared with the experimental 
graphene flakes. In the experiment, if we take the trench width as the length of graphene, 
the maximum thermal conductivity for a 10.6 nm long graphene is 5000/2000*10.6= 26.5 
W/mK, which is in the same order of the present work but still quite different. This difference 



23] Since our main focus is to study 



is probably results from the rough quantum correction, 
the isotopic doping effect on thermal conductivity, we accept this quantum correction for 
all calculations in present paper. More accurate quantum effect can be considered by other 
more precise method, such as the nonequilibrium Green's function method (see Ref. |24| for 
review) . 

Fig. [2] shows the isotopic doping effect on the thermal conductivity. The carbon isotope 
14 C is added by randomly substituting 12 C in pure 12 C GNR system, and the value of the 
thermal conductivity in the doped GNR is obtained by averaging over ten randomly doping 
processes. It shows that the thermal conductivity decreases rapidly before 10% doping, and 
decreases slowly after this doping percentage. With only 10% isotopic doping, the therma 
conductivity is reduced up to 50%, which is comparable with the results in nanotubes. [a. Il4| 

In order to understand the underlying mechanism, it is essential to study localized phonon 
modes around the isotopic doping center. We apply the GULP to obtain the full force con- 
stant matrix of the GNR under the same boundary condition as used in the MD simulation. 
Then the force constant matrix is diagonalized to obtain the full phonon dispersion and 
eigen vector of the system. The localization property of each phonon mode can be analyzed 



from its corresponding eigen vector. If the component of one atom in the eigen vector is 
very large and decreases exponentially to its neighboring atoms, then this is the localized 
phonon mode around this particular atom. 

Fig. [3] shows the phonon modes localized around one 14 C doping center. There are three 
localized modes around this 14 C doping atom, with frequencies as 634.4, 1373.3 and 1629.9 
cm -1 . Due to their localizing property, these localized modes have negative effect on the 
thermal conductivity. The frequencies of the first two localized mode are not very high. 
They can be excited in room temperature region. So these localized modes can reduce the 
thermal conductance effectively. 

When two 14 C doping centers occur together in the GNR, there will be two different 
cases. 

(1). In the first case, these two 14 C doping atoms are far from each other. Fig. H] shows 
that the frequencies of these three localized modes do not change. Since these two doping 
centers are far from each other, they can not 'see' each other and these modes can be 
localized around either of these two 14 C atoms, leading to degeneracy of all three modes. 
This degeneracy means that there are two phonon modes which have the same frequency. 
However, from the eigen vector of them we can see that they localized around different 
doping atoms. So in this case the number of localized modes will increase quickly (linearly) 
with increasing doping atoms. 

(2). In the second case, these two 14 C atoms are very close with each other. Fig. [5] 
indicates that these two 14 C atoms act more like a single C-C molecule doping center. 
Their effect is to slightly change the frequency of the first two modes, but keep them non- 
degenerate. While the third mode becomes degenerate without changing frequency. So in 
this second case the number of localized modes increases very slowly with increasing doping 
atoms. We stress that the number of the first two low frequency localized modes does not 
increase. These two modes are more important than the third high frequency localized mode 
to reduce thermal conductivity. 

We have also studied localized modes for more than two isotopic doping centers. The 
situation is similar as two doping centers described above. In the first case where isotopic 
doping centers are far from each other, the number of localized modes increases linearly. 
While in the second case, this number increases slowly. 

We can now explain the isotopic doping effect on the thermal conductivity in Fig. [2J In 



the low doping region (< 10%), most isotopic doping atoms are far from each other, which 
is the first case mentioned above. So the number of localized modes increases rapidly with 
increasing doping. As a result, the thermal conductivity decreases quickly. In the high 
doping region most 14 C atoms are close with each other. This is in the second case and 
the number of localized modes increases slowly with further increasing doping. Especially, 
the number of low frequency localized modes, which are more important in the reduction of 
thermal condctivity, does not increase too much. Consequently, the decreasing behavior of 
thermal conductivity with increasing doping slows down. 

The temperature dependence for the thermal conductivity of pure and 20% doped GNR 
is shown in Fig. El The common phenomenon in pure and doped GNR is that the thermal 
conductivity increases with increasing temperature. There are two reasons for this phe- 
nomenon. (1). At higher temperature, more high frequency phonon modes will be excited, 
which is beneficial for thermal transport. (2). As shown in Fig. (TJ the quantum correction 
is more serious in lower temperature, which reduces the value of the thermal conductivity. 
The difference between these two system is that the thermal conductivity in the doped GNR 
increases much more slowly with increasing temperature than that in the pure GNR. This is 
because localized modes in the doped GNR are more excited at higher temperature. Thus 
slow down the increasing behavior of thermal conductivity with increasing temperature. We 
stress that the increase of thermal conductivity is only true in the graphene nanoribbon 
where the edge effect is important. In large flakes of graphene, the Umklapp scattering 
dominates the phonon transport, and leads to decreasing thermal conductivity with the 
increase of temperature. 

III. CONCLUSION 

In conclusion, we study the thermal conductivity of GNR by means of molecular dy- 
namics with quantum correction. It shows that the isotopic doping can reduce thermal 
conductivity remarkably by localizing phonon modes. Since the number of localized modes 
increase linearly with increasing doping in the low doping region, the reduction of thermal 
conductivity is distinct in the low isotopic doping region. But in the high doping region, the 
number of localized modes increase slowly with increasing doping, thus thermal conductivity 
decreases slowly with further increasing doping. It also shows that the thermal conductivity 
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will increase with increasing temperature and isotopic doping can slow down this behavior 
by localizing phonon modes. 
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FIG. 2: Thermal conductivity with quantum correction v.s isotope percentage at 300 K for GNR 
with length 106 A and width 4.92 A. Line is a guide to the eye. 



E 

o 

i— 
-Q 
> 



0.6 
0.4 
0.2 







(0=634.4 



(a) 



^Wv~ 



(0=1373.3 



(b) 



(0=1629.9 



(c) 



0.2 0.4 0.6 0.8 

x/L 



FIG. 3: (Color online) Normalized vibration amplitudes vs. reduced x for three phonon modes 
localized around one isotopic doping center in GNR with length 106 A and width 4.92 A. The 
frequency u for each mode is in cm -1 . 
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FIG. 4: (Color online) Normalized vibration amplitudes vs. reduced x of three localized modes in 
GNR with two isotopic doping centers, which are far from each other. All modes are degenerate. 
The frequency ui for each mode is in cm -1 . 
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FIG. 5: (Color online) Normalized vibration amplitudes vs. reduced x of three edge modes in 
GNR with two isotopic doping centers, which are close with each other. Modes in (a) and (b) are 
non-degenerate, while the mode (c) is degenerate. The frequency u for each mode is in cm" 1 . 
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FIG. 6: (Color online) The dependence of thermal conductivity on temperature for pure and zu ■< 
doped GNR. 
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